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The discovery of the jet quenching in central Au + Au collisions at the Relativistic Heavy- 
ion Collider (RHIC) at Brookhavcn National Laboratory has provided clear evidence for 
the formation of strongly interacting dense matter. It has been predicted to occur due to 
the energy loss of high energy partons that propagate through the quark gluon plasma. 
In this paper we investigate the dependence of the parton energy loss due to elastic 
scatterings in a parton plasma on the value of the strong coupling and its running with 
the evolution of the system. We analyze different prescriptions for the QCD coupling 
and calculate the energy and length dependence of the fractional energy loss. Moreover, 
the partonic quenching factor for light and heavy quarks is estimated. We found that the 
predicted enhancement of the heavy to light hadrons (D/tt) ratio is strongly dependent 
on the running of the QCD coupling constant. 

Keywords: Quark gluon plasma; Collisional energy loss; Running coupling constant 
PACS numbers: 12.38.Mh; 13.85. N; 25.75.-q 



1. Introduction 

Relativistic heavy ion collisions provide an opportunity to study the QCD proper- 
ties at ener gy densities about thirty times higher than the density of the atomic 
nuclei U-i^isa j n these extreme conditions a deconfined state of quarks and gluons, 
the Quark Gluon Plasma (QGP), is expected to be formed in the early stage of 
the collision. Recently, the discovery of jet quenching in central Au + Au collisions 
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at the Relativistic Heavy-ion Collider (RHIC) at Brookhaven National Laboratory 
has provided clear evidence for the formation of strongly interacting dense matter. 
Detailed analysis indicate that the suppression of the single hadron spectra at high 
Pt, the disappearance of back-to-back correlation of high px hadrons and the az- 
imuthal anisotropy of high px hadron spectra in noncentral collisions are caused by 
parton energy loss (For a recent review see, e.g, Ref. El). Basically, the high parton 
density produced in heavy ion collisions could induce a large amount of energy loss 
while hard partons produced in the initial stage of the collision propagate through 
the fireball, due to the interactions of the hard partons with the medium. 

The total energy loss of a particle in a medium can be decomposed into a colli- 
sional and a radiative contribution. While the first one originates from the energy 
transfer to the medium particles, the latter one is caused by radiation from the 
fast particle. At large energies one expects that radiative energy loss becomes much 
larger than the collisional one, as in the electromagnetic case. However, at lower 
energies these two processes can contribute equally, with the collisional one being 
the dominant for small values of the parton energy. Currently, an open question 
is to quantify the contribution of each process in the RHIC kinematical region. In 
particular, in the last few years the understanding of parton energy loss by gluon 
bremss trahlung has been extensively developed (For recent reviews see, e.g., Refs. 
| 5 | 6 | 7 | 8 |j j^ ecen t works have analyzed t he energy color charge and mass depen- 
dence of radiative parton energy loss I H | 12 | 13 l^ ag we ^ as non _Abelian feature 

A basic characteristic from these works is the assumption that the radiative 
energy loss dominates disre gardin g the collisional one. Howeve r , rece nt studies in 
charm quark thermalizationHSHSI ) quenching of hadron spectra EHH anc i the elas- 
tic parton energy loss including all 2 — * 2 processes indicate that in the RHIC 
kinematical region it is far from clear that radiative energy loss dominates over 
collisional energy loss. 

Another important aspect is that in general the approaches for radiative and 
collisional energy loss start from the assumption that the properties of the medium 
and its interactions with the energetic parton project ile do not change with time 
(For a discussion of the medium evolution see, e.g. Refs.^HH^ However, in nucleus - 
nucleus collisions at collider energies, the produced hard partons propagate through 
a rapidly expanding medium. The density of scattering centers is expected to reach 
a maximum value at the plasma formation time, To, and then decrease with time r 
rapidly due to the strong longitudinal expansion. We begin following the progress 
of the quark at r = tq, where it begins a process of scattering and diffusion in the 
plasma which causes a loss of the initial momentum and relaxation towards the 
thermal velocity. The scenario which we assume in this pa per for the equilibration 
has been proposed some years ago ^2! (g e e also Ref. ^3) and consider that the 
interactions in the plasma of quarks and gluons get stronger with the time because 
the average parton energy drops due to the expansion of the system. The basic 
idea comes from the feature that the QGP is not a static medium, but it is cooling 
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while partons propagate through. Consequently, the scale of the coupling changes 
with the evolution of the system, and this feature motivates the calculation of the 
observables assuming the running coupling. In this paper we will assume that the 
temperature is the dominant scale and consequently will control the running of 
the QCD coupling. The standard procedure is to take the solution for the running 
coupling from the renormalization group equation, which to lowest (1 - loop) order 
is given by: 

MS 

and put the renormalization point ~p proportional to the first Matsubara frequency, 
~p oc 27rT, which for massless quarks is the only dimensionfull quantity inherent to 
the theory. Moreover, the strong coupling a s also depends on the scale parameter 
of the modified minimal subtraction scheme Aqcd- However, it is not evident that 
temperature is the relevant scale at the energies currently accessible in heavy ion 
physics, where the temperatures reached are only moderately larger than the phase 
transition temperature. As the temperature drops over the lifetime of the QGP, a s 
should also vary during the equilibration and the evolution of the plasma. In order to 
simplify our considerations we model the space-time evolution of the quark-gluon 
plasma by the Bjorken scenario wi th boost invariant longitudinal expansion and 
conserved entropy per rapidity unit^H This implies that T = To(^ L )3, where To is 
the initial temperature at the initial time To. We neglect the transverse expansion 
of the system. Furthermore, we will consider in what follows two prescriptions for 
the temperature dependence of the running coupling constant. First we will assume 

asiT > = (33-2n / )ln[(19T c /A 17 ^)(T/T c )]' (2) 

where T c / Aj^g — 1.78±0.03. This prescription, which we denote by thermal [ai" 1 ^], 
has been used in, e.g., Refs. I ^5 | 26 | Furthermore, we also consider that a s can be 
given by 

fT s 21)95 

a *( T ) = 7 7T~\ ~ f / 7, \T ( 3 ) 



2 In 



At 



with Q = 2ttT . This parameterization of the strong coupling has been obtained from 
recent results in the lattice ^4 We denote this prescription by a^ at ^ . 

In Fig. [T] we present the temperature dependence of a s predicted by these two 
prescriptions. We use T — 375 MeV and to = 0.33 fm which are reliable values for 
RHIC energies^^. The horizontal lines characterize the constant values of the QCD 
coupling: (a) a s — 0.3, as in previous calculations^^, (b) a s — 0.2, which is a typical 
value for the coupling constant and (c) a s = 0.5 for comparison. We can see that 
there is a large difference in normalization and shape between the two prescriptions. 
Comparing the results, the lattice motivated prescription given in Eq. ([3]) [ai 1 " 1 ^ in 
the figure] provides the largest a s value at small temperatures. The fixed values 
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a s = 0.2 and 0.3 lie between these two prescriptions. On the other hand, a fixed 
value of a s = 0.5 is larger than the lattice prescription at large tem peratu res. It is 
important to emphasize that if the prescriptions proposed in Refs. | 22 | 28 | are usec j 
we obtain a behavior for the running coupling similar to ai lat ^. 

As the temporal de velopment of th e coupling constant modifies the behavior of 
several QGP signatures ^ l ^^^^ , we can expect a similar effect in the estimates 
of the energy loss of a parton propagating in a QGP. Our main goal in this paper is 
to estimate the influence of the running coupling cons tant in the collisional energy 
loss, presenting a reanalyzes of the studies from Refs. | 17 | 18 | 31 |32|^ anc j an estimate 
of the quenching factor Q{pi_) for light and heavy quarks (See discussion in Sect. 
[2]). We postpone for a future publication the study of this effect in the radiative 
energy loss (For a short discussion see Ref. &^). 

A comment is in order here. In our studies we consider that the fracti onal en ergy 
loss is given in terms of the Bjorken formula ^31 as generalized in Refs. I 35 | 36 [ j^ e _ 
cently, Peshier^^has advocated that the collisional energy loss is an observable for 
which loop corrections to the tree level approximation are essential. In particular, if 
the running of the coupling constant is taken into account the collisional energy loss 
becomes independent of the jet energy in the high energy limit {E 3> T). It implies 
a distinct temperature dependence for the elastic energy loss when compared to 
the Bjorken formula and a larger value for the mean energy loss per length. If the 
formalism considered in this paper is reanalyzed assuming as input the Peshier's pre- 
diction, we can expect a larger contribution of the elastic energy loss. Consequently, 
the results presented here can be considered as a lower bound for the contribution 
of the elastic energy loss for the parton quenching. 

This paper is organized as follow. In next section we define the partonic quench- 
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Fig. 1. Temperature dependence of the QCD coupling considering two different prescriptions. 
Sec text. 
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ing factor Q(p±) and present its relation with the energy loss AE. Furthermore, the 
Fokker-Planck equation is derived and its solution is obtained. In Sect.[3]we present 
our predictions for the time evolution of the drag coefficient for the propagation 
of light and heavy quarks considering different prescriptions for the QCD coupling, 
as well as our results for the energy and length dependence of the fractional en- 
ergy loss. Moreover, our predictions for the quenching factor for light and heavy 
quarks are also presented. Finally, in Sect. [4] we present some remarks and our main 
conclusions. 



2. Partem Propagation in a QGP and the Partonic Quenching 
Factor 

Lets start our discussion considering the hadron production in nucleus-nucleus col- 
lisions. Due to the interactions of the produced parton with the medium, it loss an 
additional energy fraction AE while escaping the collision. Consequently, the in- 
clusive transverse momentum spectra of the particles produced in nucleus - nucleus 
collisions will be modified with respect to hadron - hadron collisions. In general, 
the hadron formation is described in terms of parton recombination and/or by the 
fragmentation of the energetic partons. In particular, it is expected that for p± > 5 
GeV the hadrons are dominantly produced by fragmentation. In what follows we 
assume that fragmentation is the dominant process of hadron formation and that 
it occurs after the parton has left the comoving medium. It allows to analyze the 
energy loss effects directly in the p± spectrum of the scattered partons. In Refs. 
.38 .39 the e ff ec t f the radiative parton energy loss in the p± distribution has been 
estima ted. Here we extend this approach for collisional energy loss. Following Refs. 
| 38 | 39 | we assume that the p± spectrum is given by: 

~¥pY = J deD(e) — rfV — = Qip±) *p± ' (4) 

where dN ji^^ +e ^ is the transverse momentum distribution in elementary parton- 
parton collisions, evaluated at a shifted value p± + e, and D(e) is the probability 
distribution in the energy e lost by the partons in the medium through collisions. 
Moreover, Q(p±) is the medium dependent quenching factor. In Ref. 1^81 the authors 
have demonstrated that in a realistic calculation of the quenching, the knowledge 
of the full probability distribution is actually required. However, as our goal is 
to study the effect of the running of coupling c onst ant in the estimates of the 

on 

quenching factor, we will assume, following Ref. , that the quenching can be 
modeled in terms of the mean energy loss. Therefore, in this work we wi ll ca lculate 
the partonic quenching factor assuming that it can be approximated by ^21 

dN mcd dN vac 
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where 

dN mcd _ 1 f 2 '^ f R 2 ^ dN vac (p ± + AE) 



d*p ± 2w*B* J ^ J dr d 2 P± ' (6) 

R is the nuclear radius, <f> is the angle between the velocity and the radius vector of 
the parton [See Eq. (fT5|) ] and AE is the total energy loss by partons in the medium. 

I n o rder to calculate AE we will use the approach proposed by Svetitsky in 
Ref. ^H, which considers the Brownian motion of a parton in a thermal bath, gov- 
erned by the Fokker-Planck equation. This approach has be en used to estimate the 
diffusion of char m qua rks in a quark - gluon plasma I 31 | 32 | ^g ee a j g0 j^ e f [43]) anc j 
its equ ilibration I 44 | 45 |^ ag we jj ag ^ calculate the quenching of light and heavy 
quarks D-5 | 16 | 18 [ j n ^his approach one starts from the Boltzmann equation for the 
distribution function f(x,p) and assume that there is no external force acting on 
the quark and that the phase space distribution / does not depend on the posi- 
tion of the quark. Furthermore, assuming the scattering process to be dominated 
by small momentum transfers, one arrives at a Fokker-Planck equation describing 
the evolution of / in the momentum space. We then investigate the time evolution 
of the Fokker-Planck equation in a thermally evolving QGP. The main simplifying 
assumption in the approach is the momentum-independence of the drag and dif- 
fusion coefficients. In what follows we present the main formulae of this approach. 
However, we do not repeat here t he details of the calculations and refer the reader 
to the original works IM1 32 | 17 | 18 | 

The Boltzmann equation in the relativistically covariant form can be written as 

M 

p»d fi f(x >P )=C{f}, (7) 

where p^ = (Ep } p) is the four-momentum of the test quark and / is its phase-space 
density and C{/} is the collision term. Following Refs. |3U 32 | 17 | 18 | we assume that: 
(a) the hydrodynamical evolution can be described by the Bjorken scenario, which 
implies that it is valid to assume that the plasma is uniform and consequently the 
phase space density of the quark is independent of x, and (b) the collision term is 
given by the elastic collisions of the test quark with other quarks, antiquarks and 
gluons in the system. Using the Landau app roximation and restringing the analyzes 
to the one-dimensional problem one obtain I 31 | 32 | 17 H8| 

df d d 2 

m = dp- [Tlip)f] + w 

which is the Landau kinetic equation, with the transport coefficients given by 
Ti(p) = J dkw(p,k)k=^ = (F), 

Up) = \ Jd 3 kw(p,k)k 2 = { -^, (9) 

where w(p, k) sums the rate of collisions of a test particle with the partons of the 
medium and (F) is the average force acting on the test particle. If we consider 



i = zPi(p)fl + n^. (8) 
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that the background heat bath is constituted of a large amount of weakly coupled 
particles in thermal equilibrium at a temperature T, with some non-thermal but 
homogeneously distributed particles due to the fluctuations, the problem can be 
simplified SSI assuming that the equilibrium of the bath will not be disturbed by 
the presence of these few non-thermal particles. Due to its small number, one can 
also assume that they will not interact among themselves, only with particles of the 
thermal bath. Consequently, one can replace the distribution functions of the colli- 
sion partners of the test particle by their Fermi-Dirac or Bose- Einstein distributions 
and Eq. ([5]) reduces to the Fokkcr-Planck equation. 

The transport coefficients 71 and T2 can be expressed in terms of the drag 
coeffi cient o f quark A, which is almost independent of momentum p as shown in 
Refs. 1211221 it allows to write the Fokker-Planck equation as follows l21l22llxIlSl 



df 
dt 



V, 



£f 
dp 2 



(10) 



where T>p is the diffusion coefficient, which is given by T>p = AT 2 if we assume 
that the momentum p can be approximated by the temperature T of the system and 
the coupling between the Brownian particle and the bath is weak ESI The Eq. (fT0|) 
describes the evolution of the momentum distribution of a test particle unde rgoing 
Brownian motion. The solution from Eq. (fit)]) has been obtained in Refs. I 17 | 1<B [ 
assuming as b oun dary conditio n / (p, t ) — —> 5(p — pn) and using the method of 
characteristics S3. It is given by | 17 | 18 | 



exp 



Pae 



■ I L W) dt' 



W(L) 



(11) 



f(p,L) = 

where W(L) is given by 
W(i) = ^4jT V F (t')exp 

' ' (12) 

and is the probability distribution in momentum space. In the previous equation 
the length of the expanding plasma is assumed as being the maximum time limit. 

From Eq. (|11|) we can estimate the mean energy of the parton due to elastic 
collisions after traversing a distance L. It is given by 









2 f A(t")dt" 




exp ^-2 J A(t') dA 


Jo 







Ef(p,L)dp 



(13) 



Then the average energy loss due to elastic collisions in the medium will be given 

by 



AE = E n 



:e\ 



(14) 



where E = m± = \/p\ + M 2 at the central rapidity region, y = 0. Consequently, 
in order to estimate the average energy loss it is necessary to calculate Eq. (jlip in 
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terms of the drag coefficient of a quark. Following Refs. we approximate the 

drag coefficient by its average value, 



,MpM; = .41/. = <-~ 



(15) 



which is directly dependent of the energy loss rate dE/dL. The abo ve approximation 
is reasonable up to moderate momentum values (p < 15 GeV) H^l. 

The energy loss rate in the QGP due to elastic collisions with high-momentum 

1041 

transfer ha ve been originally estimated by Bjorken <^ and recalculated in Refs. 
| 40 | 35 | 36 |41j taking m to account the loss with low-momentum transfer dominated by 
the interactions with plasma c ollec tive modes in the hard thermal loop approxima- 
tion S3. In particular, in Ref. ^^Sl the authors have estimated the energy loss for 
heavy quarks and in Ref. ^51 for light partons. For heavy quarks and in the domain 
E « M 2 /T, it reads 



dE 

1l 



8na 2 T 2 



1 1 



2v 2 



In 





In 









2 6+ "/ B(v) 



ET 

TTlgM 



whereas for E » M 2 /T, it is 



dE 
dL 



8na 2 T 2 



6 



hi 



2 3(6+ ™^0.92- 



(16) 



(17) 



where n/ is the number of quark flavors, a s is the strong coupling constant, m g = 
+ rif/6)g 2 T 2 /3 is the thermal gluon mass, E is the energy and M is the mass of 
the quark. B(v) is a smooth velocity function, which can be taken approximately as 



q 7I36) p or ijght quarks we use the expression (fT7|) and set M = in the calculations. 
Inclusion of the diagrams other than the t channel increase the energy loss rate for 
light quarks by a factor 2 0^2. 



3. Results and Discussion 

At the energies (temperatures) which we are interested in this paper the drag coeffi- 
cient A for partons propagating in a plasma can be calculated using Eq. (fT5|) and the 
expressions for elastic energy loss given by Eqs. (fl6|) and (fT7|) . The average over the 
momentum is made using the Boltzmann distribution. As we assume the Bjorken 
scenario for the hydrodynamics evolution ^21, the time dependence of the temper- 
ature is given by T(t) = i^To/i 1 / 3 , where t and T are, respectively, the initial 
time and temperature at which the background of the partonic system has attained 
local kinetic equilibrium. The time dependence from the drag coefficient is directly 
associated with this evolution for the temperature, which decreases with time as the 
system expands. We assume as maximum time limit for the evoluti on the length of 
the plasma L. These approximations has been considered in Refs. which we 

would like compare our results. Moreover, we assume Tn = 375 MeV and in = 0.33 
fm for RHIC energies as in Ref. The results for the time dependence of the drag 
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Fig. 2. Time evolution of drag coefficient for light (left panel) and heavy (charm) quarks (right 
panel) . 



coefficient A are show in Fig. [5] for light and heavy (charm) quarks. We present a 
comparison among the results for fixed a s and running coupling from lattice and 
thermal QCD prescriptions. One can see that the drag coefficient is very sensitive 
to the modifications in the prescription used for a s . We have that when the lattice 
one is considered (dotted line), the coefficient becomes larger than earlier results 
for fixed a s = 0.3 (solid line) in all the evolution of the fireball. In comparison to 
the a s = 0.5 case (dashed line), the lattice prescription implies a smaller value of 
A(t) for small values of t. At large times, oi£ becomes larger than 0.5 due to the 
cooling of the system, which implies a larger value of A(t) . On the other hand, when 
the thermal QCD prescription is considered (dot-dashed line), the drag coefficient 
is smaller than the previous result, being closer to the prediction obtained assuming 
a s = 0.2 (long-dashed line). 

In Fig [3] we present the results for the fractional energy loss, as a function of 
the distance traveled by the partons (i.e, the time evolution of the plasma) for light 
and heavy quarks. Following the results for the drag coefficient, shown in Fig. [2] and 
discussed above, we have that a s = 0.5 implies the larger amount of energy loss at 
small values of L. On the other hand, the lattice QCD a s prescription implies a large 
energy loss for large L. The thermal QCD a s prescription provides the smallest one. 
Moreover, as expected, the fractional energy loss for heavy quarks is smaller than 
for light quarks. 

In Fig. [4] we present the ratio between the fractional energy loss for heavy and 
light quarks. While for fixed a s the ratio is almost constant in the range considered, 
for running a s with the lattice prescription, the ratio is monotonously increasing, 
and the value is greater than 0.5 in all range considered; on the other hand, using 
the thermal prescription, the ratio is strongly suppressed when compared with the 



February 2, 2008 3:16 WSPC/INSTRUCTION FILE mackedanz'coll 



0.2 ■ 
0.1 - 



~i — 1 — i — 1 — i — 1 — r 

T s 375 MeV 
t„ = 0.33 fm 






1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 






T = 375 MeV 


a =0.2 






t = 0.33 fin 


a =0.3 






M = 1.5 GeV 


4 = 0.5 
























a , ( ' 






























































/ »* 















5 

L (fm) 



Fig. 3. Fractional energy loss as a function of the distance traveled by the light (left panel) and 
heavy quark (right panel). 
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Fig. 4. Heavy-to-light ratio of the fractional energy loss. 



result for fixed a s . This feature suggests that the heavy quarks lose less than 20% 
of the energy lost by light quarks in its path through the fireball. The predictions 
obtained using a s = 0.5 are similar to the lattice one. 

In Fig. [5] we present the results for the fractional heavy quark energy loss as 
a function of the incident quark energy. One have that the predictions are almost 
energy independent in the energy range E ~ 5 — 15 GeV for all prescriptions ana- 
lyzed. This featu re is d ue to the momentum independence of the drag coefficient, as 

lOI IQOI 

discussed earlier 11 '. Furthermore, the magnitude of the energy loss is strongly 
dependent on the prescription used in the calculation. 

In order to compute the p^-spectra, we assume that the geometry of the heavy 
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ion collision is described by a cylinder of radius R and the parton moves in the 
transverse plane in the local rest frame. Consequently, if we consider a parton cre- 
ated at a point r with an angle (f> in the transverse direction it will travel a distance 

m 

L(4>) = {R 2 - r 2 sin 2 0) 1/2 - r cos^ , (18) 
where cos0 = v • f ; v is the velocity of the parton and r — |r|. Furthermore, we 
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Fig. 5. Heavy quark fractional energy loss as a function of the energy of incident quark. 
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Fig. 6. Quenching factor for light quark px _s P ec trum. An estimative of radiative energy loss is 
presented for comparison. The left panel presents results with v = 8.0 in the vacuum spectrum 
parametrization, while results with v = 12.42 are presented in the right panel. 
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assume the following parameterization of the p± distribution 



dNjr = A ( i \ 



d 2 p± \Po+P±) 



(19) 



wher e tw o sets of parameters are available in the literature: v — 8.0 and po = 1.75 
GevESand v = 12.42 and p Q = 1.71 GeV 48 . The results for the quenching factor 
for light quarks are shown in the Fig. [5] for both sets of parameters. For comparison, 
we present the estimate of the quenchi ng due to the radiative energy loss, following 
the parametrization proposed in Ref. Due to the smaller drag coefficient, the 
thermal QCD a s prescription gives a high quenching factor, so the spectrum is less 
modified by collisional energy loss than in the case of fixed a s . At a s — 0.3 and 
v = 8.0 we have that in the high p± region, elastic and radiative energy loss are of 
the same order of magnitude. On the other hand, the lattice QCD a s prescription 
implies a very large modification of the spectrum, similar to the predictions obtained 
assuming a s = 0.5. Finally, if the system presents a lower value of a s than considered 
in earlier calculations, the gluon bremsstrahlung becomes the dominant mechanism 
of energy loss again. 

For heavy quarks we use the p± distribution of cha rme d hadrons (D-mesons) 
produced in hadron collisions. It is experimentally found to be well described by 
the following simple parametrization 



where b — 1.4 ± 0.3, n = 10.0 ± 1.2 and M c = 1.5 GeV. The quenching factor for 
heavy quarks is shown in Fig. [Jj The results are similar to those obtained for light 
quarks in Fig. [6] While the thermal QCD a s prescription gives the higher quenching 
factor, the a s = 0.5 one gives the smaller factor. The lattice QCD prescription 
implies a small quenching factor Qh ~ 0.1, shown a strong suppression in the 
charm spectra when this prescription is used. Again, at smaller values of a s , higher 
the quenching factor due to elastic scattering in the QGP. 

Recently, the ratio between the spectra of hadrons with heavy quarks and with 
light quarks has been proposed as a tool to investigate the medium formed in heavy 

[To] 

ion collisions Because of its large mass, radiative energy loss for heavy quarks 
would be lower than for light quarks. It occurs due to combined mass effects HHU2J 
the reduction of the formation time of gluon radiation and the suppression of gluon 
radiation at angles smaller than the ratio of the quark mass to its energy by destruc- 
tive quantum interference^^ - the dead-cone effect. The predicted consequence of 
these distinct radiative energy losses is an enhancement of the heavy to light spectra 
ratio at moderately large transverse momentum, relative to that observ ed in the ab- 
sence of energy loss (A recent analysis for LHC energies is given in Ref. . As the 
behavior of this ratio considering collisional energy loss is still an open question, in 
Fig. [8] we present the ratio between heavy and light quark quenching factors, which 
reflects the heavy to light hadrons (D/tt) ratio, considering only collisional energy 
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loss. While for values of fixed a s smaller than 0.5, the results show an enhancement 
factor close to 1.4, this enhancement is suppressed when running coupling prescrip- 
tions are considered. This feature could suggest that with an expanding cooling 
medium, the collisional energy loss for heavy and light quarks would be of similar 
magnitude. Finally, if a large value of a s is considered the ratio can be smaller than 
one. 



4. Conclusions 

Before we summarize, let us discuss the assumptions made in this work in order to 
simplify the calculations, which may affect our results. First, we have disregarded 
the momentum dependence of the drag and diffusion coefficients, which contain 
the dynamics of the elastic collisions, replacing it by its average value. Second, the 
entire discussion is based on the one dimensional Fokker-Planck equation and the 
Bjorken model for a nuclear collision, which may not be a very realistic description, 
providing only a qualitative estimate. An extension for three dimensional analysis 
is still in discussion and may lead to a revision of our conclusions. Moreover, the 
inclusion of diagrams for other channels in elastic processes implies an additional 
uncertainty. The perturbative expressions for the (radiative and collisional) partonic 
energy loss, which are only known to leading order in the strong coupling constant, 
may get substantial corrections of higher order. A calculation of next-to-leading 
order corrections to energy loss in perturbative QCD would be desirable to test the 
stability of our results and conclusions. Finally, in order to compare our predictions 
with the experimental data the fragmentation of partons into hadrons should be 
considered and included in our calculations. 

As a summary, in this paper we have investigated the dependence of the parton 
collisional energy loss in a QGP on the value of the strong coupling. Since the plasma 
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Fig. 8. Heavy-to-light quenching factor ratio. 

is not a static medium, a fixed value for a s is a crude approximation. More realistic 
estimates should take into account the evolution of the fireball. We have considered 
running coupling in the calculation, evolving it with the cooling of the QGP. From 
the Fokker-Planck equation, we derived the transport coefficients and related them 
with parton mean energy loss. The drag coefficient is found to be modified with 
the value of a s considered, and it is strongly dependent on the running coupling 
prescription used in the analysis. The lattice one gives larger drag coefficient values, 
while the thermal QCD one gives smaller value, always compared with the fixed 
value used in previous calculations. A similar feature is verified in the mean energy 
loss results and quenching factors. For light quarks, we found that the radiative 
and collisional energy loss are of the same order of magnitude, in the high p± 
region, if the a s value is larger than 0.3. For smaller values of coupling, the gluon 
bremsstrahlung becomes the dominant process for energy loss, again. The ratio 
between heavy and light quenching factors has been studied, and we have found 
the absence of enhancement if running coupling is used. It was a striking result, 
since it suggest that heavy and light quarks have the same order of magnitude for 
the collisional energy loss. Our results motivate a similar study in radiative parton 
energy loss. 
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